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—_ Grcces d=) Abstract: In this study, we proposed some inflation 
rate predictions based on econometric models that 

performed better than the targets of the National Bank 
of Romania. Few econometric models (multiple 
regressions model and a vector-autoregression) were 
used to predict the quarterly inflation rate in Romania 
during 2000:Q1-2016:Q4. The GDP growth has a 
negative impact on inflation rate in Romania, an 
increase in logarithm of GDP with one percentage 
point determining a decrease in inflation logarithm 
with less than 0.1 units according to both proposed 
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pe models. However, an increase in inflation rate in the 
previous period determined an increase in this 
variable in the current period. The inverse of 
DOI: unemployment rate is positively correlated with the 
10.14254/jems.2016.1-1.1 index of prices. The causal relationship between 


inflation rate and unemployment rate is reciprocal. In 
the first period the index of prices evolution is 
explained only by changes in this variable. The 
inflation rate volatility is due mainly to the evolution 
of this indicator, the influence decreasing 
insignificantly in time, not descending under 88%. 
More than 99% of the variation in unemployment rate 
is explained by the own volatility for all lags. The 
annual forecasts based on these models performed 
better than the targets on the horizon 2015-2016. 
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unemployment rate. 


1. Introduction 


An important objective of the Romanian monetary authorities is to maintain the inflation 
rate in a certain target. Before providing forecasts, it is important to know the past evolution of this 
variable. Falnita and Sipos (2007) proposed a multiple regression model for inflation rate in 
Romania in the enlarged EU.Damian (2011) studied the disinflation process in Romania within the 
context of the European integration. The problem of forecasting the inflation rate in the new 
members of EU was also analyzed by Arratibel, Kamps and Leiner-Killinger (2009). The authors 
showed that the recent rise in inflation was due to adverse supply shocks, such as energy and food 
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price increases, as well as to ongoing changes in the economic structures of these catching-up 
economies. 

The main aim of this study is to propose better annual inflation rate than the official target. 
We construct several econometric models to explain the evolution of quarterly index of prices, 
using multiple regression models and a VAR model (vectorial-autoregressive model) The VAR 
approach allows us to evaluate the variance decomposition of each indicator. In this way, we can 
determine if the variation in the variable’s evolution is mainly due to the other variable or to its 
own evolution. Woglom (2000) also used VAR model for explaining quarterly inflation that was 
affected by external shocks in the South Africa in 90’s. The author obtained a small linkage between 
inflation rate and monetary policy. Therefore, the inflation targeting was not recommended for 
South Africa. 

We use as forecast method only simple econometric models, avoiding complex models, 
because according to Armstrong, Green and Graefe (2015), the golden rule in forecasting is to be 
conservative. Moreover, for inflation rate in G-7 countries, Canova (2007) showed that the 
univariate model performed better than the multivariate ones. 

Inflation rate is an important macroeconomic indicator used by many institutes in decision- 
making process. The central bank is directly interested by the most accurate inflation rate forecasts 
in its targeting. Abubakar (2016) made a critical review of the arguments in favour or against 
inflation targeting in order to ensure the price stability in a country. We consider that inflation 
targeting is a good tool in establishing the monetary policy, but the values of the targets should be 
more realistic. Our opinion is also consistent with previous studies in literature of Ball and 
Sheridan (2004) and Goncalves and Salles (2008). The authors showed that for emerging countries 
like Romania, inflation targeting did not succeeded in reducing the inflation volatility. Our solution 
is not to eliminate the inflation targeting, but to provide better values for targets that take into 
account the structure of the economy. 

The proposal of some targets for inflation to support the monetary policy strategy is an 
option that was taken into account also by Bernanke and Woodford (1997) who showed the 
importance to determine if actual policies are consistent with long-term inflation target. The 
authors showed that the Central Bank might use as inflation targets the predictions of private 
forecasters, but it also has to take into account the structural model of the economy in order to 
guide the policy decisions. Moreover, Stein (2015) noticed that the inflation rate expectations are 
affected by the seasonality of price index that causes deviations in measuring the inflation 
phenomenon and in predicting it. 

In this article, after the methodological approach regarding the inflation modeling, several 
empirical models are proposed for explaining the quarterly index of prices evolution in Romania 
during 2000:Q1-2014:Q4. 


2. Modeling inflation rate 


Some studies analyzed the pass-through of external shocks, like shocks in exchange rate, oil 
price, and non-oil import price to inflation at different stages of distributions. The analysis 
employed vector-autoregression processes (VAR model) that included the repartition chain of 
pricing. 

More types of models were used to explain the evolution of consumer prices index: 


ICP, - index of consumer prices at time ¢; 

IC P,—4 - index of consumer prices at time ¢; 
GCF; - real gross domestic product at time ¢; 
U; - unemployment rate at time ¢; 

€+ - error term; 

Qo, 1, Hz, Az - parameters. 


The model is linearized by the introduction of the variable U’ that is the inverse of the 
unemployment rate: 


log(ICP,) = a@) + a, logUCP,_,) + az log(GDP,) + a,U', + &, (2) 
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L'. - inverse of the unemployment rate at time t. 
Another econometric model is proposed to explain the price index: 


log(ICP,) = a0+ al * log(ICP,_,) + a3 * log(GDP,) + £1; (3) 


a0, al, a2- parameters; 

¢_1t- error term. 

We consider y_i a vector of “m” variables. Each of these variables has “p” lags. The rest of the 
variables (the deterministic variables and the constant) are placed in a vector denoted by y..; that 
has m* elements. The VAR model has the following form: 


yi = ACL) y-1 + Cy *+ ey, = AC)¥i-1 + Cy *;+ @,e; > (0,%e), (4) 
Number of regressors: k= mq+m* 
Number of coefficients: c=mk 


The VAR model is written in two equivalent forms (X - a Tk matrix, Y and E- Tm matrices, 
Im- identity matrix, a- mk vector, y and e- mT vectors ): 


Y=XA+E 
y= Um* X)at+e,e > (0, * Ir) (5) 
3. Econometric models for quarterly inflation 


We start from the relationship between gross domestic product, inflation rate and 
unemployment rate. Quarterly data were collected for Romania using the data base of National 
Institute of Statistic and Eurostat, covering the period 2000:Q1-2014:Q4. The data are seasonally 
adjusted using Tramo/Seats method. In order to ensure the data stationary, some transformations 
are made to the data. The logarithm was applied to index of consumer prices and real GDP. 

The following models were used to explain the evolution of consumer prices index: 


M1: log(ICP,) = 0.785 + 0.33 « log(ICP,_1) — 0.07 « log(GDP,) + 0.32 *U;, + & (6) 
M2: log(ICP,) = 0.472 + 0.52 « log(/CP,_,) — 0.02 « log(GDP,) + €1; (7) 


The errors are independent and homoskedastic. Moreover, we do not have reasons to reject 
the normal distribution of the errors. The GDP growth has a negative impact on inflation rate in 
Romania, an increase in logarithm of GDP with one percentage point determining a decrease in 
inflation logarithm with less than 0.1 units according to both proposed models. However, an 
increase in inflation rate in the previous period determined an increase in this variable in the 
current period. The inverse of unemployment rate is positively correlated with the index of prices. 

A VAR model is used to explain the price index and unemployment rate evolution in Romania 
for the same period (2000:Q1-2014:Q4). 

The Granger causality test is applied for data series in order to establish if a variable is cause 
for the other one. In Granger acceptance, a variable X is cause for Y if better predictions result when 
the information provided by X is taken into account. 

The results of Granger causality test show that unemployment rate is the cause of inflation 
rate and the inflation rate is the cause of unemployment. 


Table 1: VAR Granger causality tests 
Dependent variable: LOG_CPI 


4.629431 


0.156699 . 
All 0.156699 0.9246 


Source: author's computations 


4.629431 
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Almost all the lag length criteria, excepting FPE, at 5% level indicate that a VAR(2) model is 
the best model. 





Table 2: Lag length criteria 
flag | kogh =| LRT | SO FPE | AIC | SC 
os} 118.0939 | NA | 2.45E-05 | -4.940166 | -4.861436 
5.459258 
5.552518" 
5.376724 
5.246355 


Source: author's computations 








All the tests necessary to be applied for checking the validity of the estimated VAR(2) model 
are displayed in the following tables. The form of the VAR model is the following: 


LOG_CPI = 0.2061948594*LOG_CPI(-1) + 0.5444928221*LOG_CPI(-2) - 0.00686234259*U(-1) - 
0.001448886657*U(-2) + 0.003352129141 


U = 1.539690788*LOG_CPI(-1) - 0.6364159519*LOG_CPI(-2) + 0.369408104*U(-1) - 
0.1735908506*U(-2) - 0.01830193332 


VAR Residual Portmanteau Tests are used to test the errors’ autocorrelation for both 
identified model. The assumptions of the test are formulated as: 

HO: the errors are not auto-correlated 

H1: the errors are auto-correlated 

For the lag 1 up to 12, the probabilities (Prob.) of the tests are greater than 0.05, fact that 
implies that there is not enough evidence to reject the null hypothesis (HO). So, we do not have 
enough reasons to say that the errors are auto-correlated. So, after the application of Residual 
Portmanteau Test, the conclusion is that there are not autocorrelations between errors for VAR(1) 
model. 


Table 3: Residual Portmanteau test for errors auto-correlation 


flags | Stat | Prob. | AdjQ-Stat_| Prob. | df 
cr} 2.080013 [Naw] 2.124265 | Naw | NAY 
20 | 2.218787 | NA* | 2.268007, || EO NA* | NA* 
3 | 461288 | 0.3846 | 4.337192 | 0.3623, | 
4 | 9.62976 | 0.0118 | 24.18064 | 0.0671 | 8 
6 | 23.95599_ | 0.0905 | 26.05102_ | 0.0533, | 16 
so | 31.40307_ | 0.1425 | 34.92490 | 0.0695 | 24 
9 | 3.63316 | 0.2132 | 3 7.65676 | 0.1051 | 28 


Source: author's computations 





The homoscedasticity is checked using a VAR Residual LM test for the VAR(1) model. If the 
value of LM statistic is greater than the critical value, the errors series is heteroskedastic. LM test 
shows that there is a constant variance of the errors, because of the values greater than 0.05 for the 
probability. The Residual Heteroskedasticity test is applied in two variants: with cross terms and 
without cross terms. 
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Table 4: VAR Residual Heteroskedasticity Tests: No Cross Terms (only levels and squares) 


Jointtest: | | 
Chisg Taf Pro 
4.29765 | 4 Tze S | 


Joint test: ae ee ee ee 


Chi-sq df Prob. 


27.52404 09586 [| FT 


|Dependent___|_R:squared_| _F(14,34) | Prob. __'|_Chi-sq(14) | 
frest*rest | 0.083264 | 0.220578 | 0.9979 | 4.079929 | 
jres2*res2_ | 0.225447 | 0.706878 | 0.7523 | 1404691 | 


0.244235 0.784822 0.6777 11.96749 | | 


Source: author's computations 
















The normality tests are applied under the Cholesky (Lutkepohl) orthogonalization. If the 
Jarque-Bera statistic is lower than the critical value there is not enough evidence to reject the 
normal distribution of the errors. 


Table 5: VAR Residual Normality Tests Orthogonalization: Cholesky (Lutkepohl) 
Jomt |i 230445 08912 


Te ee OO OOO 
Joint | id (OG 8983 
a a 
|Component___| Jarque-Bera | df] Prob. | 
dof cotozzg0 29 
20 TC*dC 0.282572, BZ 
Joint sit nagesa 789 


Source: author's computations 





The Residual normality test provided probabilities greater than 0.05, fact that implies that 
the errors series has a normal distribution when Cholesky (Lutkepohl) orthogonalization is applied. 
The impulse-response analysis and the decomposition of error variance are made. 
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Figure 1: The responses of each variable to own shocks or the other variable shocks 


Response to Cholesky One S.D. Innovations + 2 S.E. 


Response of LOG_CPI to LOG_CPI Response of LOG_CPI to U 











Response of U to LOG_CPI Response of U to U 











Source: author's graph 


In the first period the index of prices evolution is explained only by changes in this variable. 
The inflation rate volatility is due mainly to the evolution of this indicator, the influence decreasing 
insignificantly in time, not descending under 88%. More than 99% of the variation in 
unemployment rate is explained by the own volatility for all lags. 


Table 6: Variance decomposition of the variables 

Variance | CSCCC*|s«S@DeecompositionofLOG_CPY | 
8.500938 
10.67199 
10.83687 
(6 to781 | 8858657 11.41343 
11.55304 
8 | toes | 88.202490 1.79751 
oT 5a 88.0905 1189095 
| Variance | CT SCieecompositionofU | 
6 T 08524500 | 0448329985517 
8 08524600] 049367 | 9.85063 
oT 08524638 | 50456 | 9.84954 
10 





10 0.352466 0.151714 99.84829 


Source: author’s computations 
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Another application could be the construction of a Bayesian VAR for these two variables. 
Mester (1998) proposed a VEC model for inflation rate in Romania. The VECM model proved that 
the basic transmission mechanism runs from base money (via interest rates which affect the 
relative return on financial assets) to the exchange rate and then to prices. 


4. Inflation rate predictions and targets 


The presented econometric models are used to make quarterly inflation rate for Romania on 
the horizon 2015:Q1-2016:Q4. The quarterly predictions are aggregated by computing the mean in 
order to obtain the forecasts for the annual inflation rate. These annual forecasts are compared 
with the annual targets that are established by the National Bank of Romania. 

Table 7: Annual inflation rate forecasts, actual values and targets (%) (horizon: 2015-2017) 


Year 
Mimodel_| M2model_| VARmodel | | 





Source: author’s computations 


The inflation rate predictions and the targets for 2015 and 2016 were evaluated. The 
forecast error is computed as difference between actual value and predicted value (target). 


Table 8: Inflation rate forecasts’ errors in absolute value (percentage points) (horizon: 2015- 
VA VB Te) 


Quarter Forecasts based on: 


M1 model M2 model VAR model 
2015 
2016 


Source: author's computations 





In 2015 and 2016, the VAR model performed better than the multiple regression models in 
terms of forecast accuracy. Moreover, the targets used by Central Bank were quite far from actual 
values and from our predictions. In this context, we recommend the use of VAR model predictions 
in establishing the value of the inflation target in Romania, under the assumption that the historical 
performance will maintain also in the near future. 


5. Conclusions 


According to this analysis based on VAR model, we can conclude that for the quarterly 
inflation in Romania during 2000-2015, the causal relationship between inflation rate and 
unemployment rate is reciprocal. 

In the first period the index of prices evolution is explained only by changes in this variable. 
The inflation rate volatility is due mainly to the evolution of this indicator, the influence decreasing 
insignificantly in time, not descending under 88%. More than 99% of the variation in 
unemployment rate is explained by the own volatility for all lags. The GDP growth has a negative 
impact on inflation rate in Romania, an increase in logarithm of GDP with one percentage point 
determining a decrease in inflation logarithm with less than 0.1 units according to both proposed 
models. However, an increase in inflation rate in the previous period determined an increase in this 
variable in the current period. The inverse of unemployment rate is positively correlated with the 
index of prices. 

The proposed models were used to build quarterly inflation rate forecasts. The annual 
predictions based on these expectations were compared to the targets established by the National 
Bank of Romania. The results indicated that these predictions based on simple econometric models 
performed better than the NBR targets. Therefore, the National Bank could fix a more realistic 
target given by these forecasts. 


Appendix A. Supplementary material 


Supplementary data associated with this article can be found, in the online version, at 
http://dx.doi.org/10.14254/jems.2016.1-1.1 
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